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We use observational data from Type la Supernovae (SN), Baryon Acoustic Oscillations (BAO), 
Cosmic Microwave Background (CMB) and observational Hubble data (OHD), and the Markov 
Chain Monte Carlo (MCMC) method, to constrain the cosmological scenario of holographic dark 
energy with varying gravitational constant. We consider both flat and non-flat background geometry, 
and we present the corresponding constraints and contour-plots of the model parameters. We 
conclude that the scenario is compatible with observations. In la we find Qao = 0.721q' O 3, fli,o = 

-0.0013tg;{£to. c = °- 80 ±o:i4 and A G = G'/G = -0.0025±Ho5°> while for the present value of the 
dark energy equation-of-state parameter we obtain wo = — 1.041q'2o- 
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I. INTRODUCTION 



Nowadays it is strongly believed that the universe is experiencing an accelerated expansion, and this is supported 
by many cosmological observations, such as SNe la [l|, WMAP [2|, SDSS [|[ and X-ray Q. A first direction that 
could provide an explanation of this remarkable phenomenon is to introduce the concept of dark energy, with the most 
obvious theoretical candidate being the cosmological constant. However, at least in an effective level, the dynamical 
nature of dark energy can also originate from a variable cosmological "constant" or form various fields, such is a 
canonical scalar field (quintessence) [6[ , a phantom field, that is a scalar field with a negative sign of the kinetic term 
, or the combination of quintessence and phantom in a unified model named quintom Q . The second direction that 
could explain the acceleration is to modify the gravitational theory itself, such is the generalization to /(i?)-gravity 
9]), scalar-tensor theories with non-minimal coupling [lOj . string- inspired models fllj etc. 

Going beyond the aforementioned effective description requires a deeper understanding of the underlying theory 
of quantum gravity unknown at present. However, physicists can still make some attempts to probe the nature of 
dark energy according to some basic quantum gravitational principles. Currently, an interesting such an attempt is 
the so-called "Holographic Dark Energy" proposal [12|, ll3j . Its framework is the black hole thermodynamics and the 
connection (known from AdS/CFT correspondence) of the UV cut-of of a quantum field theory, which gives rise to 
the vacuum energy, with the largest distance of the theory [l4j |. Thus, determining an appropriate quantity L to serve 
as an IR cut-off, imposing the constraint that the total vacuum energy in the corresponding maximum volume must 
not be greater than the mass of a black hole of the same size, and saturating the inequality, one identifies the acquired 
vacuum energy as holographic dark energy: 

3c 2 

PA = 8^GP' (1) 

with G the Newton's gravitational constant and c a constant. The holographic dark energy scenario has been tested 
and constrained by various astronomical observations 15] and it has been extended to various frameworks 

EMU- 

However, there are indications that Newton's "constant" G can by varying, being a function of time or equivale ntly 
of the scale factor In particular, observations of Hulse- Taylor binary pulsar [20, [21J, helio-seismological data [221 . 
Type la supernova observations [l| and astereoseismological data from the pulsating white dwarf star G117-B15A [291 ] 
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lead to G/G ^ 4.10 x 10 u yr 1 , for z < 3.5 [2J|. Thus, in our previous paper [241 ], we investigated the holographic 

dark energy scenario under a varying gravitational constant, and we extracted the corresponding corrections to the 
dark energy equation-of-state parameter. 

In the present work we are interested in constraining the model of holographic dark energy with a varying Newton's 
constant, using various observational data arising from SN, BAO, CMB and OHD. Such an analysis is crucial for 
the validity of the aforementioned scenario. The plan of the work is as follows: In section [TT| we present briefly the 
holographic dark energy scenario with a varying Newton's constant and the corresponding expressions of the dark- 
energy equation-of-state parameter. In section HTT1 we perform a combined observational constraint analysis, allowing 
for variations in all model parameters. Finally, section ITVl is devoted to the summary of our results. 



II. HOLOGRAPHIC DARK ENERGY WITH VARYING GRAVITATIONAL CONSTANT 

In this section we briefly review the holographic dark energy proposal, in the case where the gravitational constant 
is itself a function of the scale factor [241 ] ■ In order for our results to be more transparent we examine separately the 
flat and non-flat background geometries. 

A. Flat universe 



In the case where the space-time geometry is a flat Robertson- Walker: 



ds 2 = -dt 2 + a(t) 2 (dr 2 



2 dn 2 ), 



with a(t) the scale factor and t the comoving time. As usual, the first Friedmann equation reads: 



H* = 



IttG 



Pa 



(2) 



(3) 



with H the Hubble parameter, p m = ^p 2 -, where p m and p A stand respectively for matter and dark energy densities 
and the index marks the present value of a quantity. Furthermore, we will use the density parameter Oa ~~ s "' 
which, imposing explicitly the holographic nature of dark energy according to relation ([I]), becomes 



3H- 



PA, 



n A 



H 2 L 2 ' 



(4) 



As usual, in the case of a flat universe, the best choice for the definition of L is to identify it with the future event 
horizon [ll, G]| S3 , that is L = R h {a) with 



R h (a) = a 



df 



da' 



lt o(f) J a Ha' 2 
In this case, using ^ and the Friedmann equation ([3]), one can show that 

2 % /0 A " 



n' A = n A {i-n A ) 



i 



-0 A (l-fi A )A G , 



(•5) 



(6) 



where primes denote the derivatives with respect to In a. In this expression the first term is the usual holographic 
dark energy differential equation [l3j], while the second term is the correction arising from the varying nature of G, 
which is quantified by the parameter = G' /G. Finally, concerning the dark-energy equation-of-state parameter 
u>, considered as a function of the redshift z, we use the approximation [26| 



w(z) W Wo + Wl 



l + z 



with 124 



w 



Wl 




iv^a - <w ( i + -v^) - (1 '";* a g . 

be \ c / be 



(7) 

(8) 
(9) 
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In these expressions, the index stands for the value of a quantity at present, where the Hubble parameter is Hq and 
the scale factor is a = 1. 

In summary, using the results of this subsection, and the usual relations Q m = Q m oa~ 3 (with f2 m = |^p m the 
matter density parameter) and f2 A = QaoO>~ + > we can write the Friedmann equation ([3]) in a form suitable for 
observational elaboration as: 



H 2 = Hi 



{n m o(i + ^) 3 + o AO (i + ^ 1+ ^} 



(10) 



B. Non-flat universe 



Let us now generalize the above result in the case of a general FRW universe with line element 

dr 2 



ds 2 = -dt 2 + a 2 (t) 



1 — kr 2 



r 2 dn 2 



(11) 



in comoving coordinates (t,r,0,tp), where k denotes the spacial curvature with k = —1,0, 1 corresponding to open, 
flat and closed universe respectively. The first Friedmann equation writes: 



H z + — 

a 1 



8nG 



(prn + p\j- 



In this case, the cosmological length L in @ is considered to be jl6j |: 

L = 

where 

1 



Sinn 7^ — , 

V «(*) J 



(12) 



(13) 



sinn(V|fc|y) 



sin y k = +1, 

y k = 0, 

sinhy k = — 1. 



(14) 



A straightforward calculation using (UJ) and (fl"2")l leads to [24] : 



1 — fifc — fl\ H cosn 



n A (i - Ok - n A ) 



G_ 
G ' 



(15) 



where 



cosy k = +1, 
cosn(s/ \k\y) = {l k = 0, 

coshy k = — 1. 



(16) 



In expression (| 15|) we have also introduced the curvature density parameter fi^ = — r^jpp ■ We mention that since in 
this work we focus on observational constraints, we have adopted the minus-sign convention for fi^-definition, which 
is the usual one in observational works. Thus, in (fT5"|) J7 fc has the opposite sign comparing to 24|. Clearly, for k = 
(and thus fij,. = 0) it results to ([5]). Finally, concerning the dark-energy equation-of-state parameter w(z) we acquire 

w(z) kwq + wx ( ifi) , with [H| 



w 



1 2 

3 ~ 3c 



— + ^V^A + c 2 n k0 \i - n k0 - n A0 + - (1 - Ovo) \/n A0 + c 2 o fc0 j - ^v&c^(i - n 

3 DC L C J DC 



(17) 

fco -«Ao)A G . (18) 



In summary, using the results of this subsection, we can write the Friedmann equation (|12[) in a form suitable for 
observational elaboration as: 



H 2 = H, 



l {f! fe0 (i + z? + 0*0(1 + zf + n A0 (i + z f 1+w ^} 



(19) 
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III. OBSERVATIONAL CONSTRAINTS 



Let us now proceed to a combined observational constraint analysis of the scenario at hand, allowing for variations 
in all the aforementioned model parameters. We use observational data from Type la Supernovae (SN), Baryon 
Acoustic Oscillations (BAO), Cosmic Microwave Background [Z^(z*), R(z*), z*] (CMB) and Observational Hubble 
Data (OHD). The precise methods are summarized in the Appendix. 

In our calculations we take the total likelihood L oc e~ x I" 1 to be the product of the separate likelihoods of SN, 
BAO, CMB and OHD. Thus, the x 2 is 

X 2 {Ps)=X 2 SN+X 2 BAO+XcMB + XoHDi ( 20 ) 

and the parameter vector reads 

Ps = {n b h 2 ,n c h 2 ,n k ,c,A G }. (21) 

In addition, we obtain the three derived parameters SI AO: ^mo — ^6 + an d the Hubble constant Hq, based on the 
above basic cosmological parameters. In our analysis, we perform a global fitting on determining the cosmological 
parameters using the Markov Chain Monte Carlo (MCMC) method. The MCMC method is based on the publicly 
available CosmoMC package [27j . which has been modified to include the codes about BAO, CMB [Z^(z*), R(z*), z*] 
and OHD. Finally, apart from these and the two independent model parameters c and Aq, the basic cosmological 
parameters are taken in the following priors: the present physical baryon density Qbh 2 £ [0.005,0.9], the present 
physical cold dark matter energy density Sl c /z 2 £ [0.01,0.99], and for the non-flat case, the additional parameter 

n k £ [-0.1,0.1]. 

Using these techniques we are able to impose constraints on the various parameters of the scenario of holographic 
dark energy with varying gravitational constant. 



A. Flat universe 



For the case of flat background geometry the cosmologically interesting parameters are the dimensionless quantities 
S1ao 5 c and Aq = G'/G, while we have to take into account the additional uncertainty in Hq. The corresponding 
ID and 2D likelihood-contours are depicted in Fig. [T] We mention that for completeness we provide these plots not 
only for the parameters used in the fits, as discussed above, but also for the derived ones, although some of them are 
related to each other. 

Additionally, in order to provide the results in a more transparent way, in Table U we present the la best- fit values of 
the, used in the fits and derived ones, model parameters. 



Xmin 


flAo(lc) 


c(lcr) 






477.5 


n 79 o+0.026 
U. 1 ^O_ Q 3 Q 




-0 0016+ 0049 
U.UU-LO_ 0049 


70.21*; 



TABLE I: The minimum value of \ 2 an d the lcr best-fit values of the model parameters, in the flat-background case. All 
quantities are dimensionless, apart from Hq which is measured in km s _1 Mpc -1 . 



These figures show that the scenario at hand can be compatible with observations. Furthermore, note that Aq = 
G'/G is restricted around zero, in a region which is in agreement with independent observations and estimation of 
G'/G. In particular, observations of Hulse- Taylor binary pulsar B1913 + 16 lead to the estimation G/G ~ 2±4 x 
iQ-i2y r -i [2(|[2l| i while helio-seismological data provide the bound —1.6 x 10~ 12 yr _1 < G/G < [22^ . Similarly, 
Type la supernova observations give the best upper bound of the variation of G as — 10~ 11 yr^ 1 < § < [2Sj . while 
astereoseismological data from the pulsating white dwarf star G117-B15A lead to ^ < 4.10 x 10 _11 i/7 ,_1 [29(. (See 

also [23| for various bounds on G/G.) Since the limits of G-variation are given for G/G in units yr _1 , and since 
G/G = HG'/G, we can estimate their implied Ac substituting the value of H in yr~ l . Thus, inserting an average 
estimation for the Hubble parameter H « (H) w 6 x lO^yr" 1 [H, we obtain that < |A G | < 0.08 [II]. In 
summary, as we observe from Fig. [TJ A G at 2a is well inside these bounds, and this offers a self-consistency test for 
our analysis. 

For completeness, we close this subsection by estimating the la bounds of the dark energy equation-of-state pa- 
rameter. In particular, since we have approximated it as w(z) ps wq + uq I j^— ) , with u> and «q given by (JS]),©, we 
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FIG. 1: (Color online) The ID and 2D likelihood plots, in the flat-background case, using SN, BAO, CMB and OHD observa- 
tional data. The curves stand for the la and 2a regions. All quantities are dimensionless, apart from Ho which is measured in 
km s" 1 Mpc -1 . 



can easily calculate that in la: 

<n =i). I V :;■;'„ (22) 



w = -1.04+ - 13 



0.12 
0.04 



Interestingly enough, we observe that in this scenario the best-fit value of the present value of w, namely wq, is 
smaller than the corresponding one of simple holographic dark energy model [l3j , and it lies in the phantom regime. 
Furthermore, note that the left and right bounds are larger as expected, since we have the extra freedom in varying 
the gravitational constant. Finally, W\ is slightly larger comparing to simple holographic dark energy scenario. 



B. Non-flat universe 



We use the combination of observational data from SNIa, BAO and CMB to construct the likelihood contours for 
the free model parameters, which in this case are the dimensionless quantities Qao, c, Aq and ttko, together with Ho. 
In Fig. [2] we present the corresponding ID and 2D likelihood-contours, not only for the parameters used in the fits, 
but also for the derived ones. Additionally, in Table [TT1 we present the la best-fit values of the model parameters. 
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FIG. 2: (Color online) The ID and 2D likelihood plots, in the non-flat scenario, using SN, BAO, CMB and OHD observational 
data. The curves stand for the la and la regions. All quantities are dimensionless, apart from Ho which is measured in 
km s" 1 Mpc" 1 . 



Xmin 


SI ao {la) 


c{la) 




tt k0 {la) 


H>{la) 


477.4 


n 79 +0.03 
u - ' z -0.03 


o.8olS:3 


-n nn2 c i +0 ' 0080 

u.uuzo_ 0050 


-n nm , H +0 ' 0130 

u.uuio_ 0040 


70A± 3 2 ° 9 



TABLE II: The minimum value of \ 2 an d the la best-fit values of the model parameters, in the non-flat scenario. All quantities 
are dimensionless, apart from Ho which is measured in km s -1 Mpc^ 1 . 



As we observe, the scenario of holographic dark energy with varying gravitational constant in a non-flat background 
geometry, is compatible with observations. As expected, the extra freedom in varying the gravitational constant leads 
to larger bounds of all the parameters, comparing to those of usual holographic dark energy models [lf|. However, as 
we see from the above analysis, the variation range for Qko is relatively narrow comparing to the physically interesting 
range arising from the usual investigations of the literature [3l| . This feature has to be considered as a disadvantage 
of the present scenario. 
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Since we have approximated the dark energy equation-of-state parameter as w(z) m wq + w\ f j^J , with wo and 
wi given by (TTT)) . (fT5)l . we can easily estimate that in la: 

wo = -lMt° fo 

wi = 0.15tHo- (23) 

As we see, the present value of w, namely Wq, can lie in the phantom regime. Additionally, note that the best-fit 
values are very close to those of the flat scenario, since the best-fit value for ilko is very close to zero. However, the 
left and right bounds are larger as expected. Finally, note that the best-fit values of both wo and wi are in very good 
agreement with the corresponding ones of usual (constant-G) holographic dark energy models [l5j], but the bounds 
are larger due to the extra freedom in the G-variation. 



IV. CONCLUSIONS 



In this work we confronted the scenario of holographic dark energy with varying gravitational constant, with data 
from SNIa, CMB, BAO and OHD observations. We performed independent fittings for flat and non-flat background 
geometries, and we extracted the corresponding observational constraints on the free model parameters of each case. 

In the case of a flat universe, we deduced that the scenario at hand is compatible with observations, and we presented 
the corresponding contour-plots on the free parameters, namely f^Ao, c and Ac, taking into account the uncertainty 
in Ho- The gravitational constant variation is inside the corresponding bounds acquired from independent estimations 
of the literature, and this acts as a self-consistency test of the model. The corresponding lex-bounds are presented 
in Table HI Moreover, the best-fit value of the present dark energy equation-of-state parameter is smaller than the 
corresponding one of simple holographic dark energy models [l5j], lying inside the phantom regime (w = — 1.042+^124)- 
Finally, the left and right bounds are larger as expected, since the scenario at hand possesses the additional freedom 
of varying the gravitational constant. 

For the non-flat geometry, we showed that the scenario at hand is compatible with observations. We constructed 
the corresponding likelihood contours of the free parameters fl\o , flko , c and Aq , taking into account the uncertainty 
in Ho, and for clarity we presented the lcr-bounds in Table HT1 Furthermore, the best-fit value of the present dark 
energy equation-of-state parameter lies inside the phantom regime (wq = —1.042^0 204)3 an( i ^ is in agreement with 
the corresponding one of usual holographic dark energy models [151 ]. However, as expected, the additional variation 
of the gravitational constant leads to an increase on the extreme values. 

In summary, we conclude that the scenario of holographic dark energy with varying gravitational constant can be 
a candidate for the description of dark energy. 
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Appendix: Observational data and constraints 



In this appendix we briefly review the main sources of observational constraints used in this work, namely Type 
la Supernovae, Baryon Acoustic Oscillations, Cosmic Microwave Background [Ia(z*),R(z*),z*] and Observational 
Hubble Data (OHD). 



1. Type la Supernovae constraints 



We use the 397 SN la Constitution dataset, which includes 397 SN la [32|]. Following 
corresponding constraints by fitting the distance modulus fj,(z) as 



one can obtain the 



15 G 
fHh (*) = 5 log 10 [D L (z)] + — log 



10 



eff 
G 



Mo, 



(A.l) 
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where G is the current value of effective Newton's constant G e //. In this expression Dl(z) is the Hubble- free 
luminosity distance Hgd^z) / c, with H the Hubble constant, defined through the re-normalized quantity h as H = 
100/i km s _1 Mpc _1 , and 

/i = 42.38 - 51og 10 ft. 



where sinnn(^ jOfcJx) respectively denotes sin(y|f2fejx), ^/|f2fc|a;, sinh(-y/|f2fc|x) for Sl^ < 0, rZ^ = and f2/- > 0. 
Additionally, the observed distance moduli fJ, bs{zi) of SN la at Zj is 



Vobs(zi) = m b s {z t ) - M, 



(A.2) 



where M is their absolute magnitudes. 

For the SN la dataset, the best-fit values of the parameters p s can be determined by a likelihood analysis, based 
on the calculation of 



x 2 (p w) = ^° hs ^ zi " > " ^(poi^y 



SN 



07 



E 

SN 



{5log w [D L (p s ,Zi)] - m ob s{zi) + M'} 2 



(A.3) 



where M' = fj,o + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically i:35ij as 



resulting to 



with 



X 2 ( Ps ) = -21n 



exp 



B 



X = A - — + In — 



-\x 2 {Ps,M') 



C 



dM', 



C 



2?r 



(A.4) 



{5 \og 10 [D L (p s ,Zi)} - m obs (zi)} 2 



SN 



^og 10 [D L (p s ,Zi)] -m obs {zi) 



SN 



SN 



Relation (IA.3I) has a minimum at the nuisance parameter value M' = B/C, which contains information of the values 
of h and M. Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, note 
that the expression 

X 2 SN (Ps)=A-(B 2 /C), 



which coincides to (|A.4[) up to a constant, is often used in the likelihood analysis [33|, [35|, [36| , and thus in this case 
the results will not be affected by a flat M' distribution. 



2. Baryon Acoustic Oscillation constraints 



The Baryon Acoustic Oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, and measure the distance-redshift relation at z = 0.2. Additionally, Baryon Acoustic Oscillations in the 
clustering of the SDSS luminous red galaxies measure the distance-redshift relation at z = 0.35. The observed scale 
of the BAO calculated from these samples, as well as from the combined sample, are jointly analyzed using estimates 
of the correlated errors to constrain the form of the distance measure Dy(z) [37H39| 



D v {z) 



{l + zfD\(z) 



H(z) 



1/3 



(A.5) 
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In this expression Da{z) is the proper (not comoving) angular diameter distance, which has the following relation 
with di,{z) 

D A ( Z )- dL[z) 



(1 + z) 2 

The peak positions of the BAO depend on the ratio of Dy(z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) zj, which can be obtained by using a fitting formula [40| : 

1291 (Q h 2 )~o.4i9 

+ (A ' 6) 

with 

b 1 = 0.313(r> m /i 2 )-°- 419 [l + 0.607(r> m ^ 2 ) - 674 ] (A.7) 

b 2 = 0.238{n m h 2 ) 0223 . (A.8) 

In this work we use the data of r s {zd)/ 'Dy(z) extracted from the Sloan Digital Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) [39], which are listed in Table [TOl with r s (z) the comoving sound 
horizon size 

( ) f* c s dt _ c f a c s da _ c f°° dz ^^ f 1/il+z) da 

C Jo a C J a 2 H C J Z Z H(z) ^3 J a 2 H(a) y/1 + (3Q 6 /(4Q 7 )a) ' 

where c s is the sound speed of the photon— baryon fluid [4ll - |43T | : 



3 P-y(Z) o S2 7 



and fi 7 = 2.469 x lQ- 5 h~ 2 for T C mb = 2.75A 



z 


r s (z d )/D v (z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE III: The observational r s (z d )/D v (z) data 

Using the data of BAO in Table IPTII the inverse covariance matrix V^ 1 reads (38j : 

1 = ( 30124.1 -17226.9 \ (A 
^ -17226.9 86976.6 J' [ ' 1 

Thus, finally, the Xbao(Ps) ^ s gi ven as 

X 2 BAo(Ps)=X t V- 1 X, (A.12) 
where A is a column vector formed from the values of theory minus the corresponding observational data, with 

/ r B (z d ) _ _ 1905 33 \ 



x =[M,-^y ,A ' i3) 



(A* denotes the transpose). 
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3. Cosmic Microwave Background constraints 

The CMB shift parameter R is provided by [II] 



R{z*) = 



sum 



dz' 
o H&) 



(A.14) 



where the redshift z % (the decoupling epoch of photons) is obtained using the fitting function [45 

z, = 1048 [l + 0.00124(fVi 2 )- - 738 ] [1 + gi(n m h 2 ) 92 ] , 
and where the functions g\ and g^ read 

9l = 0.0783(^ 2 )^ - 238 [l + 39.5(^ 2 )°- 763 ] _1 
ff2 = 0.560 [1 + 21. l^t/i 2 ) 1 - 81 ]" 1 . 

Additionally, the acoustic scale is related to the first distance ratio and is expressed as 



I 



r s {z*) 



:Smn 



dz' 



o H(z') 



(A.15) 



Using the data of I a, R, z„ from [3l| which are listed in Table HVt and their covariance matrix of [Ia(z*), R(z*), z*] 
presented in |3l| : 



cr 1 



f 1.800 27.968 -1.103 
27.968 5667.577 -92.263 
y -1.103 -92.263 2.923 



we can finally calculate the likelihood L as Xcmb = ~ 21nL: 

X 2 c MB (Ps) = Ad l [C- 1 (d l ,d J )}[^! 
where Ad, = di — df ata is a row vector, and di = (I a, R,z*). 



i 

»J > 



(A.16) 



(A.17) 





5-year maximum likelihood 


error, a 




302.10 


0.86 




1.710 


0.019 




1090.04 


0.93 



TABLE IV: The values of l A (z,), R{z„), and z„, from 5-year WMAP results [31 



4. Observational Hubble Data constraints 

The observational Hubble data are based on differential ages of the galaxies (46|. In [47|, Jimenez et al. obtained 
an independent estimate for the Hubble parameter using the method developed in (46[, and used it to constrain the 
equation of state of dark energy. The Hubble parameter, depending on the differential ages as a function of the 
redshift z, can be written as 

*m-iti§- (A - 18) 

Therefore, once dz/dt is known, H(z) is directly obtained [48j. By using the differential ages of passively-evolving 
galaxies from the Gemini Deep Deep Survey (GDDS) [4{| and archival data [13 , Simon et al. obtained H{z) in the 
range of < z < 1.8 [HI]. The twelve observational Hubble data from [5l| are listed in Table W\ 
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2 


0.1 0.17 0.27 0.4 0.48 


0.88 


0.9 


1.30 


1.43 


1.53 1.75 


H{z) (km s^Mpc" 1 ) 


74.2 69 83 77 95 97 


90 


117 


168 


177 


140 202 


lcr uncertainty 


±3.6 ±12 ±8 ±14 ±17 ±60 


±40 


±23 


±17 


±18 


±14 ±40 


TABLE V: The observational H(z) 


data 











Furthermore, in [52J the authors used the BAO scale as a standard ruler in the radial direction, and they obtained 
three more additional data: H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H(z = 0.43) = 86.45 ± 3.27. 
The best-fit values of the model parameters from observational Hubble data (48j are determined by minimizing 

2 / \ [H th (p s ;zj) -H obs {z t )] 2 

XHub KPs) = 2-^ 2T~\ ' ( A - 19 ) 

i=i G {Zl) 

where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, H t s 
is the observed value, a(zi) is the standard deviation measurement uncertainty, and the summation runs over the 15 
observational Hubble data points at redshifts Zj. 
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